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1. Introduction 

Multiqubit systems are the building blocks for the architecture of quantum information 
storing and processing [1]. Their efficiency depends strongly on the possibility to 
retain quantum coherence. But these quantum systems are physically embedded and 
constantly interacting with their local quantum and classical environments which affect 
coherence [2]. Thus being fundamentally open quantum systems, they dissipate and 
decohere with time constants that depend on their parameters and coupling constants 
and on the strength of perturbation by external agents. There have been many efforts 
to increase the decoherence time, either passively by increasing isolation or engineering 
the existence of decoherence free subspaces p] or actively by affecting the dynamics 
with external controls [H O [6] . 

As it has turned out, a counterintuitive and unexpected external perturbation 
with a positive role is the action by a noise source, quantum or classical. Normally 
one expects that noise would only increase decoherence in a monotonous way, namely 
increased destruction of coherence with higher levels of noise. But it has been found 
that noise may be used to isolate a quantum system, the noise-shielding effect (N.S.), 
and may influence positively or negatively the decoherence time in a monotonous or 
non-monotonous way. These are the quantum Zeno effects and the stochastic resonance 
(S.R.) or stochastic anti-resonance (S.AR.) effects [TJ El El [IDl [HI [121 Il3] . But it appears 
that the manifestation of these effects depends, apart from the detailed dynamical setup, 
on the way the system has been prepared. 

From our preliminary investigation of noise effects on two-qubit systems [13J, it 
became evident that certain classes of initial states are more sensitive than others. 
Considering this not to be accidental, we extended our study of the state preparation 
dependence of the noisy perturbations, to 3 and 4-qubit Heisenberg XY chains. We 
have found a clear correlation between the way a quantum system of qubits reacts to 
an external classical noise and the geometry of the initial states. As a tool we used a 
Master equation for modeling the dissipative and decoherence processes, and applied 
an external classical gaussian white noise as a stochastic control. Our investigation is 
based on the decoherence properties of a 2-qubit subsystem. In the case of two qubits 
we have only the bosonic environment and the external classical noise. In the cases of 
3 and 4-qubit systems we consider the extra qubits, after they are traced out, as a form 
of local fermionic environment. We quantify the entanglement with the concurrence of 
the 2-qubit subsystem. We study numerically the dependence of disentanglement time 
(or entanglement sudden death (E.S.D.)) [H [151 [lEl [HI [13 [19l [20] on the strength of 
the applied noise for zero and nonzero temperature and for various initial states. We 
have observed the effects of increase of decoherence, stochastic resonance, stochastic 
anti-resonance and noise shielding. 

The main result is that these effects depend strongly on the initial preparation of 
the compound system and the placement of classical noise. The importance of this result 
could be appreciated in the cases where the behavior of a studied subsystem depends 
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on the preparation of a bigger system whose parts have been traced out, but they are 
in interaction with the subsystem of interest. Our results are summarized as follows: 

• 2 qubits 

(i) Stochastic anti-resonance is observed if the initial density matrix contains 
population elements P22 or P33. This becomes evident in Bell states |\Ef), 
(figure El (a) (upper curve)), while the entanglement of Bell states 1$) decays 
monotonically over increasing noise strength (figure [2] (b)). 

(ii) For vanishing dissipation rate (7 — )■ 0), the Hilbert subspace which exhibits 
stochastic anti-resonance tends to become a decoherence free subspace [T3] . 

(iii) There are no noise shields in the 2-qubit case. Noise shields require the exten- 
sion of the system to more qubits, in order to apply noise on the local fermionic 
environment. 

(iv) Temperature degrades stochastic anti-resonance very quickly. An average exci- 
tation number (n) ~ 0.5 is enough to eliminate the appearance of the anti- 
resonance (figure [2] (a) (lower curve)). 

• 3 qubits 

(i) Multiple resonances (stochastic resonance and stochastic anti-resonance) are 
observed in two product states when noise affects the traced out qubit (figure [3] 
(a)). Moreover, there are two product states that exhibit N.S. behavior , and 
two that exhibit S.AR. (lATj) . 

(ii) Product states with parallel spins do not present interesting behavior. Noise 
decreases monotonically their entanglement evolution. This holds for any ar- 
rangement of noise perturbation ( JAll) . 
(iii) By altering the initial 3-qubit preparation of a given 2-qubit state, we observe 
different behaviour in the 2-qubit disentanglement time. This becomes evident 
with the reduced 2-qubit state of a W state. A small change of initial 3-qubit 
preparation results in N.S. , while W state preparation results in S.AR. (figure H] 
(a,b)). 
(iv) Bell state |$) preparations exhibit N.S. when noise affects the traced out qubit, 
for all of the initial 3-qubit preparations (figure [3] (b) and lA2l) . 

(v) Most of the Bell state |\E') preparations exhibit N.S. when noise affects the 
traced out qubit, but there is an initial system preparation which results in 
S.R. (figure El^c) andlMD. 

• 4 qubits 

(i) Product states with parallel spins exhibit the same behavior as in the 3-qubit 

case fiBTl) . 
(ii) Most of the product states exhibit noise shield when noise is environmental 

(figure [5] (a) and IB1|) . 
(iii) There are four product states which exhibit different kinds of effects, depending 
on the placement of noise and one that exhibits S.AR. behavior when noise is 
environmental ( JBll) . 
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(iv) The noise shields tend to become weaker if we increase the anisotropy of the 
system (figure [5] (b)). 

(v) Most of the Bell state |$^), |^^) preparations continue to exhibit N.S. when 
noise is environmental, except from two I'l'"'") preparations that result in 
stochastic anti-resonance, two preparations of \'^^) which result in S.A.R., one 
|\1'"^) preparation that exhibits multiple resonances and one \^^) preparation 
that exhibits S.R. (the last one when noise is internal), (figure |6l^a) and IB2[ 
lB3ll . 

(vi) The noise shields become stronger by applying noise in both of the traced 
out qubits rather than one of them (figure [6] (b,c)). Temperature degrades 
their magnitude, but does not change their shape. There is increase of to- 
lerance against temperature compared with the 2-qubit case, because of the 
environmental qubits. We can observe a clear N.S. behavior with excitation 
values up to (n) ~ 6 (figure [Tl^a)). 
(vii) Even though it is difficult to determine the aforementioned initial preparations 
which exhibit S.AR, we can control this behaviour by increasing temperature. 
Above a critical value of the latter, S.AR. disappears and we observe N.S.. 
This attributes an interesting positive role to the temperature, for the control 
of entanglement evolution (figure [Tl^b)). 

A complete list of the effects is presented in the tables of Appendix A and B. In 
the final paragraph we classify our results in terms of the observed effects and make 
comments for the interplay between the geometry of the initial system preparation and 
the placement of noise, which seems to be crucial for their appearance. 

2. XY Heisenberg model 

The general form of a N-spin XY chain (for spin 1/2 particles) with nearest-neighbor 
interaction is 

N 

HxY = 2_^ i^JxS^S^j^i + JyS^S^j^ij , (1) 

n=l 

with 5*^ = |o",^(i = x,y,z) the spin-1/2 operators, a^^ the corresponding Pauli's ope- 
rators, h= 1 and Sn+1 = Si (periodic boundary condition). The chain is ferromagnetic 
if Jj < 0, and anti- ferromagnetic if J, > 0. Here we study this chain in the presence of 
a constant external magnetic field Uq along the z-axis. Consequently, the unperturbed 
system Hamiltonian is: 

N 

Ho = HxY + J2uJoS:^, (2) 

n=l 

for iV = 3, 4. 

The XY model has been studied for many decades because of its interesting and 
unusual features [21] ■ It is an example of integrable system. While the isotropic 
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Heisenberg chain is solved with the use of Bethe Ansatz, the XY model is solved by 
means of Jordan- Wigner transformation [22] which was introduced in 1928 and applied 
to it by Lieb et al [23j in 1961. For more information see [21]. Furthermore, due to 
its mathematical simplicity, it is suitable for our purpose and provides a basic unit for 
many quantum information implementations, such as N.M.R. quantum computation 
and quantum teleportation [251 126] . 



3. Markovian Master Equation 

We employ a Markovian form of the Master equation, where the classical gaussian white 
noise ^{t) affects the magnetic field and is appropriately incorporated in the double 
commutator [27J. For a thermal bath with T > the Master equation is 

d ^ 

n=l 

N 

+ 7(n) Y. (/^[5+]p,) - M,[V;, [V;,p,]], (3) 

n=l 

where Vz is the spin operator affected by noise. It can be any of S^ or sum of these. M^ 
is the corresponding noise parameter. For example M4 means that V4 = S"!, while M34 
means that V34 = 8^ + S^. Furthermore D[S:^]ps = S:^ps{S:^y — {{S:^yS~,ps}/2. The 
coefficient 7 is the rate of population relaxation and (n) denotes the average excitation 
quanta of the bath. It depends monotonically on the temperature and it is used to 
parameterize the latter. For T = we have (n) = and for T — )■ 00, (n) — )■ 00. The 
first term describes the unitary evolution of the system, the second and third terms the 
interaction between the system and the thermal bath, while the last one is the addition 
of an external classical gaussian white noise [27] . We assumed that each qubit has the 
same interaction with the environment, something that implies constraints on the value 
of coupling constants [131 [21]. The eigenvalues of the 3-qubit Hamiltonian are: 

(4) 

double eigenvalues, (6) 



J + 


UJ± 


ViJ-^^y 


+ 3A2 




LU± 


2 




J- 


V(J + 2a;)2 


+ 3A2 



2 
while the eigenvalues of the 4-qubit Hamiltonian are: 



± A/ J2 + A2 + 2^2 +JA^ + (4^2 + 2 J2) A2 + ( J2 - 2^2)- (7) 



± Y J2 + A2 + 2u;2 _ ^^4 + (4^2 + 2J2) A2 + ( J2 - 2u;2)- (8) 

-J ± Vu;2 + ^2^ (9) 

J±Vuj^ + A2 (10) 
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±u; double eigenvalues 
fourfold eigenvalue 



(11) 
(12) 



where J 
with u - 



^^^ and A 



Jx '-'■!/ 



Assuming common 7 for all of the qubits, we work 



4, J = 0.2, A = 0.1, so that the interaction between them has not altered 
energy level separations by more than 10% from that of the non-interacting system. 
Moreover, we choose 7 = 0.01 for the weak coupling approximation {uo ^ 7) to remain 
valid. 





(a) 3 qubits with a 2-qubit subsystem 



(b) 4 qubits with two 2-qubit subsystem 



Figure 1. The two different arrangements of the unperturbed system: 2 qubit 
subsystem of a 3-qubit chain (left) and a 4-qubit chain (right) inside a common bosonic 
heat bath. The subsystem (green circles) interacts with the bosonic bath and the local 
fermionic environment (yellow circles). The main difference between them, is that in 
the case of 3 qubits all of them interact directly, while in the case of 4 qubits the 
interaction between 1 — 3 and 2 — 4 is indirect. 



4. Entanglement Evolution 

Entanglement is a fundamental property of quantum systems and lies in the heart of 
quantum information and computation theory [HEH]. One of the most difficult problems 
in contemporary information theory is to construct measures for the quantification 
of quantum entanglement. While the entanglement of a bipartite system is well 
understood, many measures have been presented for multipartite systems and most of 
them are related to pure states. This makes even more difficult the study of decohering 
systems, where initially pure states tend to become mixed, because of the constant 
interaction with their environment. That is why we choose to work with the bipartite 
entanglement, which is quantified by concurrence and gives us the advantage to describe 
both pure and mixed states [291 133 EH |32] • 

For a system described by the density matrix p, the concurrence C is 



C 



max 



%,0), 



(13) 
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where Ai,A2,A3,A4 are the eigenvalues of spin flipped density matrix R (with Ai the 
largest one), the definition of which is: 

R = p{oy®ay)p*{ay®ay). (14) 

C lies in the range [0,1]. C = corresponds to a product state, while C = 1 to 
a maximally entangled state. All of the other states inside this range are partially 
entangled. 

The basis states for our system are formed by the tensor product of S^ eigenstates 
{\^)n,\9)n}- In the 3-qubit case we have 

{|eee), \eeg), \ege), \egg), \gee), \geg), \gge), \ggg)} (15) 

while in 4-qubit case 

{|eeee), \eeeg), \eege), \eegg), \egee), \egeg), \egge), \eggg), (16) 

\geee), \geeg), \gege), \gegg), \ggee), \ggeg), \ggge), \gggg)}. (17) 

A simplifying feature of our Hamiltonian is that the dynamics of system splits 
into two independently evolving sets of equations, resulting in smaller density matrices 
and consequently smaller differential systems. For three qubits case we get the density 
matrix lA.ll and for four qubits the density matrix IB. II correspondingly. They both have 
the half the number of elements of the original density matrix and result in the following 
reduced density matrix: 



Pl,l Pi,4 

P2,2 P2.3 



,'" 



P3,2 P3,3 



T T 

Paa Paa 



, where blank entries are equal to 0. (18) 



The mixed state defined by this submatrix is called X state (non-zero elements along the 
diagonal and the anti-diagonal) and arises in many physical situations [HH [34l [35]. It 
has the advantage that includes Bell states I^E'"'^), |$^), product states and mixed states. 
The concurrence of this matrix is equal to: 



C = max{0,Ci,C2}, (19) 



where 



C^i = 2(|pli|-v/pr^), (20) 



C2 = 2{\pl,\ - VpV^), (21) 

It is evident that the reduced matrix for 2 qubits will inherit elements from the full 
system matrix which is either 8x8 (3-qubits case) or 16 x 16 (4-qubits case). This means 
that different initial preparations of the compound system can lead to the same reduced 
matrix form, something expected, but with great consequences for the entanglement 
evolution, as we will show in the next sections. 
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5. 2-qubit chain 

For the sake of completeness we provide two new diagrams related to those of \[Z\ . In 
the 2-qubit chain there is no ambiguity for the initial state preparation, since we do 
not trace out any qubit (we want to study entanglement). The system decoheres in the 
presence of a bosonic bath and external classical white noise affecting the magnetic field. 
What we found is that if the initial preparation of qubits contains the diagonal matrix 
elements P22 or P33, then the entanglement evolves non-monotonically over noise. There 
is a value of noise which causes the shortest disentanglement time and has to be avoided. 
This is what we call stochastic anti-resonance (S.AR.). The effect becomes evident in 
the Bell states \^) (figure [2] (a) (upper curve)). Bell states |$) decay monotonically, 
since their evolution begins outside the critical subspace (figure [2] (b)). Regarding the 
bath temperature, it plays a negative role for the S.AR., something normally expected. 
A small increase of temperature is enough to eliminate the shape of S.AR. (figure |2] (a) 
(lower curve)). 
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t : t . t . t 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.0 0,1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 



(a) S.AR. 



(b) Monotonous Decay 



Figure 2. Entanglement sudden death time (E.S.D.) for 2 interacting qubits, in 
the presence of a bosonic environment and external classical white noise, (a) Bell 
state 1^^). The upper curve corresponds to zero temperature. There is a clear anti- 
resonance behaviour. The lower curve shows that a small increase of temperature 
erases the anti-resonance, (b) Bell state |$^) for T = 0. There is no anti-resonance 
in this case. This shows that the geometry of the initial state is important for the 
appearance of this effect. Noise is collective, (w = 1, 7 = 0.01, J = A = 0.1 ) 



6. 3-qubit chain with T = 

The reduced density matrix has the general form: 

^2 



P' 






^2 











Ei=iP 



i+6,i 



i=3 Pi,i 2^i=3 Pi,i+'2 
1^1=3 Pi+'^,i 1^1=5 Phi 





Z^i=l Pi,i+& 




l^i=7 Pi,i 



(22) 
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Obviously it inherits elements from the 3-qubit density matrix lA.ll Different prepa- 
rations of the latter result in the same state of the subsystem. When all of the elements 
of the compound density matrix that participate in the sums of the reduced density 
matrix are equal, we call the state "balanced" (see the first preparation of Bell states 
in lA2llA3l[B2llB3l) . 

Firstly we focus on product states of 3-qubits which result, obviously, in 2-qubit 
product states. In order to study the entanglement evolution of product states, we 
calculated the area included between the concurrence graph and the time axis until 
concurrence reaches zero for the first time after its production, for different values 
of noise strength. We are interested here in the first cycle of creation-decay of 
entanglement, which has the largest amount of the latter, and not for any rebirths. 
This measure is sensitive in all kinds of behavior the entanglement production-decay 
can exhibit (smooth or oscillatory) and carries information for both of the maximum 
of entanglement production and the disentanglement time. The results are presented 
in the table \K\\ There are two non- interesting cases where all qubits are aligned (|eee) 
and \ggg)) and noise degrades entanglement monotonically regardless of its placement, 
two states which present S.AR. behaviour, two states with multiple resonances (S.R. 
and S.AR.) and two with N.S. behavior. The results show clearly that all of the above 
effects occur when the noise is applied on the environmental qubit (M^ = M3). 

Then we move on to the 2-qubit Bell states | $'*'), l^"*") and try several initial 3-qubit 
preparations that result in them. Almost all of them exhibit noise shield behaviour when 
noise lies outside our subsystem except from one. Indeed, there is a preparation of |\&"'") 
which exhibits stochastic resonance when noise is environmental (1A21 IA31 figure [3] (b,c)). 

Finally we present our results for the bipartite subsystem of W state 
\eeg) + \ege) + \gee) 



Vs 



(23) 



which has been proposed as the maximum 3-qubit entangled state and has non-zero 
entanglement across any bipartition. The W state has matrix representation 



Pw 
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(24) 



and its reduced 2-qubit density matrix is 



I 




1 1 

3 3 

1 1 

3 3 





i 








(25) 
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with initial entanglement Cmiuai — 0.667. Again, we observe S.AR. only when the noisy 
qubit is environmental. The most interesting result, though, is that if we prepare the 
same 2-qubit state, but with different elements of the 3-qubit density matrix, we observe 
N.S. behavior (figure |l](a,b)). This can be done easily by setting p2,3 = P3,2 = in the W 
state preparation. The reduced density matrix ( l22l) does not contain these elements, so 
the 2-qubit state remains the same. The ambiguity in the preparation of the subsystem 
has important consequences on its evolution. 
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(a) S.R. and S.AR. 



(b) N.S. 







































■ Preparation A3 (2) 
Preparation A3 (3) 



















































































































































































































































































































































































0.0 0.1 0-2 0,3 0.4 0,5 0,6 0,7 O.S 0,9 1,0 

(c) N.S. and S.R. comparison 

Figure 3. 3 interacting qubits. E.S.D. time for several initial states, as a function of 
noise parameter, (a) We observe the multiple resonances (S.R and S.AR.) for the initial 
product state \geg) — )■ |ye), where the S.R is very sharp and occurs with a very small 
noise strength, while S.AR. is much smoother and occurs with a larger noise strength, 
(b) N.S. behavior of the balanced Bell state |$~^). (c) Two different preparations of 
Bell state l^"*") that exhibit S.R. and N.S. correspondingly. Noise is external (affects 
the traced out qubit) in all of the above cases, ((n) = 0, cjo = 4, J = 0.2, A — 0.1, 
7-0.01) 



Dependence on state preparation of noise induced effects in multiquhit systems 
7. 4-qubit chain with T = 

The reduced density matrix in this case is: 



11 



■5-^4 








Y^4 





v^8 


Y^8 

Z^i=5 Pi,i+i 








Y^8 







Z]j=l Pi+12,j 











The main difference from the 3-qubit case is that the qubits do not interact directly 
with each other. Qubit 1 does not interact directly with 3 (same for 2 and 4). We continue 
to study the entanglement between 1 and 2, because we are interested in the behavior 
of two directly interacting qubits. 

As we can see in table IBll there are many product states where the combination of 
entanglement production and disentanglement time increases monotonically over noise, 
when the latter is environmental. This is noise shield behavior (figure |5] (a)). We 
also observe three cases where multiple resonances become evident. Two of them 
{\ggeg), \ggge)) "do not obey the rule" and exhibit resonances only when the noise 
penetrates at least one qubit of the subsystem under study. 

The behavior of Bell states 1$"'') and I^E'^) is not the same any more. While in the 3- 
qubit case Bell states l^"*") exhibit N.S. for all of the initial preparations we have studied, 
here there are some initial preparations that cause the system to exhibit stochastic anti- 
resonance, similarly to the W state case of 3-qubits flB2p . Furthermore, we found that 
the N.S. becomes weaker when the anisotropy parameter A increases (figure [5] (b)). 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



(a) S.AR. 



(b) N.S. 



Figure 4. 3 interacting qubits. (a) E.S.D. time for the state resulting after tracing 
out the third qubit of the W state |g'^9)+|ggf)+l9ee) ^ ^^ ^ function of external noise 
parameter M^. We observe S.AR behavior, (b) By altering a little bit the 3-qubit 
preparation we get a completely different behavior, a N.S.. For small values of noise, 
the disentanglement times are quite smaller than those of the previous case. Noise is 
external (affects the traced out qubit) in all of the above cases, ((n) — Q, luq — 4, 
J = 0.2, A = 0.1, 7 = 0.01) 



Dependence on state preparation of noise induced effects in multiquhit systems 



12 



Bell state 1^1')''' presents a variety of effects depending on its initial preparation and 
the placement of noise. We found N.S., S.AR. and S.R. (IB 3 1 and figure [6] (a)). In figure[6] 
(b) and figure [6] (c) we present the results for the balanced Bell state 1^1'^). It becomes 
evident that the application of noise should take place out of the subsystem, in order to 
create a noise shield. Moreover, the later is amplified when both environmental qubits 
are noisy, rather than one, something intuitively expected. 

8. Non-zero temperature and special initial Hilbert subspaces 

As we presented above, there are some initial preparations both in 3 and 4-qubit chains 
that cause the subsystem to exhibit S.AR or S.R., while other result in N.S. behavior, 
which is obviously the best result we can get by adding noise to a quantum system. 

The identification of these "special" regions of initial Hilbert space is a very difficult 
task, because of the great number of real parameters needed to describe the density 
matrix of a 3-qubit or 4-qubit system. The XY interaction Hamiltonian and our Master 
equation may simplify significantly the system of differential equations we need to solve, 
but even in this case, the number of real independent parameters for 3 qubits is 31 while 
for 4 qubits is 127. 

Moreover, the choice of working on the bipartite entanglement does not simplify the 
problem: The evolution of the reduced 2-qubit density matrix has direct dependance on 
the non participating elements of the full density matrix, something expected, because 
one solves firstly the full differential system and then traces out the desired number 
of qubits. Needless to say that if we wish to study larger systems, then the problem 
acquires extremely large complexity and becomes practically unsolvable. 

There are two ways to approach this problem: 

(i) The theoretical one: We have to understand how the initial geometry of the full 



0.3 0.4 0.5 0.6 0.7 0.8 



A-0,I 
A=0.3 
A=0.« 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 



0.7 0.8 0.9 1.0 1.1 



(a) N.S. for a product state 



(b) N.S. for different A 



Figure 5. 4 interacting qubits (a) N.S. for the product state \eggg) — ?- \eg). (b) N.S. 
of the balanced |$^) against anisotropy constant A. The bigger the anisotropy, the 
weaker the noise shield, ((n) = 0, cjq = 4, J = 0.2, A = 0.1, 7 == 0.01) 



Dependence on state preparation of noise induced effects in multiquhit systems 



13 



system provides the necessary background for the beneficial (or not) synchronization 
of two coexisting timescales (dissipation and decoherence). This could lead to the 
optimized exploitation of noise addition in order to get the best possible result, 
which is obviously a noise shield. We believe that the solution of this problem will 
come from a deeper analysis of the geometry of quantum states and the multipartite 
entanglement measures. 

The practical one: Given the complexity of the problem, it would be helpful if 
we could control the extent of these areas by means of an external macroscopical 
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(a) E.S.D. for different preparations of |$) 



(b) Noise affecting the subsystem 
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(c) N.S. comparison for different place- 
ments of classical noise 

Figure 6. 4 interacting qubits. (a) E.S.D. time for two preparations of |^^), when 
noise affects the two traced out qubits. The black squares correspond to the first 
preparation (balanced state) in table IB3I and the red dots to the fifth one. It is 
clear that the initial state subspace determines the entanglement evolution: The first 
preparation exhibits N.S, while the second one exhibits S.AR. and is characterized 
by much longer disentanglement times for small values of noise, (b) E.S.D. time for 
the balanced Bell state 1^"*"), as a function of external noise parameter. The system 
response is better when noise affects only one "internal" qubit and not both of them, 
something expected, (c) We observe a substantial change of behavior when placing 
noise out of the subsystem. The response is monotonically positive (i.e. N.S.). We 
find a better behavior (almost linear), with much longer disentanglement times, when 
both of the environmental qubits are noisy. The same hold for all of the Bell state 
|$+) preparations, ((n) == 0, wq = 4, J = 0.2, A = 0.1, 7 = 0.01) 
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(a) E.S.D. time for a Bell State 



(b) S,AR, and N.S, for special preparations 



Figure 7. 4 interacting qubits (a) E.S,D time for the balanced Bell state \^'^) as a 
function of temperature parameter (n), for different values of M34, We observe the 
robustness of N,S, against temperature, although the absolute value of Iesd becomes 
smaller as n gets bigger (something expected), (b) Another preparation of \^^) (|B3I 
(5)) which exhibits S,AR,. It is clear that above a critical {n) we get N,S, instead of 
S,AR,, something that attributes a positive role to the bath temperature, (wq = 4, 
J = 0,2, A = 0,1, 7 = 0,01) 



In the present stage of our understanding of the problem of which classes of 
states are related to noise induced effects, we think that finding an external control 
seems more straightforward. One option is to study the effect of temperature on 
the entanglement evolution, since temperature has the advantage of being an easily 
controllable macroscopical variable. In our previous work fT3] we observed the negative 
role of a small increase of temperature in two qubits of XY model (figure [2] (a)). By 
moving to a larger system (3 and 4 qubits), we expected to see shielding effects, due to 
the environmental qubits. 

This was indeed our first result: The larger the fermionic environment we have, the 
bigger the tolerance of the subsystem against noise. Although the disentanglement time 
decreases as (n) gets bigger, the results are the same for all of the states which exhibit 
monotonic behavior in T = case: Noise shield for T = remains shield for T > as 
we can see in figure [7] (a). We note again that in the case of 4 qubits the shield becomes 
"stronger" when both of the external qubits are noisy. 

The most interesting result though, is the effect of temperature increase on those 
special initial preparations that result in S.AR. behavior. In the 4-qubit case there is a 
{n) critical under which S.AR. becomes evident. Above this value, the subsystem acquires 
N.S. behavior as we can see in figure [7| (b). Temperature eliminates the first-negative 
part of S.AR. and leaves only the second-positive part (ie. N.S.). This is a very welcome 
result, because it gives us the opportunity to get the best possible noise effect, the noise 
shield, whether the subsystem begins from a "S.AR. region" or not and provides a way 
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to avoid the identification of the latter: If we observe S.AR. then we must increase the 
temperature and get N.S.! 

9. Conclusions 

In our study of a system of two interacting qubits, under the influence of a bosonic 
environment and classical noise, we had observed that the appearance of effects of non- 
monotonic dependence of the disentanglement time on the noise strength, was dependent 
on the initial preparations of the system [13] . Then the natural question that arose was 
whether the initial preparations of the local environment of the 2-qubit subsystem does 
influence the dynamical evolution. As it has turned out, the local fermionic environment 
does indeed influence the way the subsystem reacts to an external classical noise. We 
have studied several initial preparations of 3-qubit and 4-qubit chains and focused on 
the bipartite entanglement of the reduced subsystem. Even though they differ by only 
one qubit, they represent two different cases of a quantum system: In the 3-qubit chain 
all qubits interact directly with each other as we see in figure [H This is not the case 
in general for Heisenberg type n-qubit chains with n > A. Our system is modelled 
by a Markovian Master equation with two Lindblad dissipators: one representing the 
common quantum bath and one representing the classical external field. While the first 
one is related to the amplitude damping channel, the second one implies a continuous 
indirect measurement of Vz- In |36] one can see that the last term of the Master equation 
points to the quantum Zeno effect [7], which has been proposed as a reliable method for 
the protection of entanglement against environmental noise [HI |9]. 

What we have observed, is that the application of external noise in combination 
with the proper initial preparation of the compound system, has clear influence on 
the entanglement evolution of the subsystem. While many initial preparations of a 
compound system result in the same subsystem state, they may lead to very different 
dynamical behaviours of the subsystem entanglement. The 3 interesting effects we 
observed, are: 

• Noise shield: Monotonous increase of disentanglement time against noise. This is 
the most likely effect, according to our summary tables, when noise is environmental 
(both in 3 and 4-qubit chains). The traced out qubit (s) play the role of a 
local fermionic environment which shields the subsystem against dissipation and 
decoherence. It is strongly related to quantum Zeno effect. 

• Stochastic anti-resonance: Non-monotonic behavior of the disentanglement time, 
firstly decreasing and then increasing. This means that a moderate value of noise 
has the worst possible effect on the system and needs to be avoided. This was 
initially observed in our 2-qubit study [13] for various initial states. In our 3 and 4 
qubit study, the effects depend not only on the preparation but also on the way the 
noise is applied. We observed S.AR. for 4 initial preparations in the 3-qubit chain 
case, with noise acting only on the traced out qubit. In the 4-qubit case, there were 
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9 appearances of S.AR., 7 with noise acting only on the environmental qubits and 
2 with noise inside the subsystem. 

• Stochastic resonance: Non-monotonic behavior of the disentanglement time, firstly 
increasing and then decreasing. There is a moderate value of noise that upgrades 
system reaction. We observed it 3 times in 3-qubit chains, with noise acting on the 
traced out qubit and 5 times in 4-qubit chains, 3 with noise inside the subsystem 
and 2 with noise out of the subsystem. 

We note that product and entangled states have different time evolution. Product 
states begin at zero entanglement, they evolve to entangled states because of qubit 
interaction and finally loose their entanglement in the presence of environmental 
decoherence-dissipation. We studied the area under the concurrence graph of the first- 
largest cycle of entanglement, in order to gain information for both the magnitude of 
entanglement production and the disentanglement time. Indeed, there is no unique 
measure of stochastic resonance-like effects and one has to decide which is the best for 
the physical property under study. 

Moreover, temperature can play a constructive role for the control of the special 
initial Hilbert subspaces which lead to non-monotonic noise effects. It provides us with 
an easy way to cut out the first-part of S.AR.. 

Of course the number of different preparations we tried is not large, but it is enough 
to show the main results, which strengthen our belief that: 

(i) Noise effects are not fully understood yet, even in the Markovian regime which 
simplifies significantly the mathematical calculations. 

(ii) Initial system preparation (and hence the local geometry of the state space) may 
allow noise to affect decoherence in a positive or negative way. 

(iii) The initial local state space of qubit systems needs to be mapped according to 
the noise effects it can exhibit. The classification of these subspaces could lead to 
predictions of non-monotonic behavior, and consequently to effective optimization 
of a quantum computation process. 

Appendix A. 3-qubit chain density matrix and summary tables 

The general form of the 3-qubit density matrix is: 



Ps 



(A.l) 
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where all blank elements are equal to zero. 

Next we present summary tables for the 3-qubit case, for several initial preparations 
of product states and Bell states. 



Table Al. Product state results 



Initial State 


N. S. 


S. 


AR. 


S. R. 


Noise 


|eee) -^ |ee)i2 


\eeg) -^ |ee)i2 


/ 








Ms 


lege) -^ leg) 12 




/ 
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\eg9) -^ leg) 12 




/ 
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Ms 


|gee) -^ |ge)i2 
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M3 


Iseg) -^ |.ge)i2 




/ 




/ 


Ms 


\9ge) -^ \gg) 12 


/ 








M3 


Iggg) -^ \gg) 12 
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Table A2. Bell state |$+) - I'-^'^yy' 










Initial State 




N. S. 


S. AR. 


S. R. Noise 


Pia = P2.2 = pTj = /08,8 = 1/4 

Pi, 7 = PT.l = P2..& = PS,2 = 1/4 




/ 




Ms 


Pia = 0, P2,2 - 1/2, P7,7 = 0, P8,8 = 1/2 
Pi, 7 = /07,1 = 0,/32,8 = P8,2 = 1/2 




/ 




Ms 


P14 = 2/5,p2,2 = 1/10, P7,7 = 2/5,P8,8 = 
Pi, 7 = P7,l = 2/5,/92,8 = P8,2 = 1/10 


a/10 


/ 




M3 


pi,l = 1/10, P2,2 = 2/5,p7,7 = 1/10, /?8,8 - 
Pl,7 = P7,l = 1/10, P2,8 = P8,2 = 2/5 


= 2/5 


/ 




M3 


pi,l = 1/2, /92, 2 = 0,/97,7 = 1/2, P8, 8 = 
Pl,7 = P7,l = 1/2,P2,8 = P8,2 = 




/ 




M3 



Table A3. Bell state |*+) == M+^ 



Initial State 



N. S. S. AR. S. R. Noise 



P3,3 = P4,4 = P5,5 = P&fi = 1/4 
P3,5 = P5,3 = P4,6 = P6,4 = 1/4 



/ 



M. 



P3,3 = 0,P4,4 = 1/2, P5, 5 = 0,P6,6 = 1/2 
P3,5 = P5,3 = 0,P4,6 = P6,4 = 1/2 

P3,3 = 1/2,P4,4 == 0,/35,5 = 1/2, /36, 6 
P3,5 = P5,3 = 1/2, P4,6 = P6,4 = 

P3,3 = 2/5, /94,4 = 1/10, P5,5 = 2/5, P6,6 = 1/10 / 
P3,5 = P5,3 = 2/5,p4,6 = P6,4 == 1/10 

P3,3 = 1/10, /?4,4 = 2/5, p5, 5 = 1/10, P6, 6 
P3,5 = P5,3 = 1/10, P4,6 = P6,4 = 2/5 



/ 



Ma 



= 



/ 



Ms 



Ms 



= 2/5 / 



Ms 
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Appendix B. 4-qubit chain density matrix and summary tables 

The general form of the 4-qubit density matrix is: 



(B.l) 



where all blank elements are equal to zero. 

Next we present summary tables for the 4-qubit case, for several initial preparations 
of product states and Bell states. 
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Table Bl. Product states 






Initial State 


N. S. 


S. AR. 


S. R. 


Noise 


eeee) — > ee)i2 


\eeeg) -^ ee)i2 


/ 






M34,M4 


\eege) -^ |ee)i2 


/ 






M34 


\eegg) -^ |ee)i2 


/ 






M34,Af4 


\egee) -> |ee)i2 


/ 






M34,M4 


\egeg) ->■ \eg)i2 




/ 




M34,M4 


\egge) ->■ \eg)i2 


/ 






M34,M4 


\eggg) -^ \eg)i2 


/ 






M34,M4 


\geee) -)• 56)12 


/ 






Af34,M4 


\geeg) ->■ \ge)i2 


/ 






M34,M4 


\gege) ->■ \ge)i2 


/ 


/ 




N.S.^ M4, S.AR. -> A'f34 


\gegg) -^ \ge)i2 


/ 


/ 


/ 


N.S -^ Af34, S.AR. + S.R. -^ Mi 


\ggee) -^ \gg)i2 


/ 






Msi, Mi 


\ggeg) -^ \gg)i2 




/ 


/ 


S.AR. + S.R. ^ M1234, 
M234, M12, Mi3, 
Mi4, Ml 



\999e) -^ \gg)i2 / 



/ / N.S ^ M4 

S.AR. + S.R. ^ M1234, 
M234, or M12, Mi3, 
Mu, Ml 



\9999) -^ 155)12 



S.AR. + S.R. means that there exist resonance and anti-resonance in the same 
diagram (i.e. figure |3]^a)). 
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Table B2. Bell state |$+) = 


\ee) + \99) 
V2 










Initial State 




N. S. 


S. AR. 


S. R. 


Noise 



Pia = P2.2 = P3,3 = PiA — 8 

Pl3,13 = Pl4,14 = PlbSb = PWM = § ^ -^34, ^i 

Pl3,l = Pl4,2 = Pl5,3 = Pl6,4 = g 

Pia3 = P2M ^ P3,15 = P4,16 = 8 

P14 = 5; /02,2 = P3,3 = P4,4 = 

Pi3,i3 ~ 2TPiiM — Pi5a5 = pieae ==0 / M34 

Pl3,l = Pi, 13 = 2 ' ^14,2 = Pl5,3 = 
^2,14 — Pz.lb — Pi.16 — Pl6,4 — 

Pl4 = 0,/92,2 = 5,/53,3 = P4,4 = 

Pi4,i4 =: |,/3i3,i3 = /O1545 = Pieae = / M34, M4 

n „ „ 1 



Pl3,l — Pi, 13 — 0,P2,14 — Pl4,2 — 2 
P3,15 = Pl5,3 — P4,16 — Pl6,4 — 



Pl,l = P2,2 = 0,/93,3 = p4,4 = \ 

Pl3,13 = Pl4,14 = 0,/9i5,i5 = pi6,16 = I: ^^ M34, M4 

Pl3,l — Pi, 13 = P2,14 = Pl4,2 = 0, 

P3,15 — Pl5,3 — P4,ie = Pl6,4 = 4 

Pl,l = P2,2 = |, P33 = P4,4 = 

Pl3,13 = Pl4,14 = |,Pl5,15 = Pl6,16 = / M34, M4 

Pl3,l = Pi, 13 = P2,14 = Pl4,2 = 4 

P3,15 = Pl5,3 ~ P4,16 = Pl6,4 = 

Pl,l = P2,2 = P3,3 = 0,P4,4 = 5 

Pl3,13 = Pl4,14 — pl5,15 = 0,/?i6,16 =2 ^ -^34 

Pl3,l = Pi, 13 = P2,14 = Pl4,2 = 

P3,15 == Pl5,3 — 0,p4_i6 — Pi6,4 = 2 

Pl,l = P2,2 = ie,P3,3 = P4,4 = Yg 

Pl3,13 = Pl4,14 ~ Yg,Pl5,15 ~ Pie.ie =16 '^ M34,M4 

Pl3,l = Pi, 13 = P2,14 ~ Pl4,2 = Yg 

3 
P3,15 = Pl5,3 ~ P4,16 = Pl6,4 = Jg 
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|eg) + |ge) 



Table B3. Bell state |*+) = ^^^ 



Initial State 



N.S. S.AR. S.R. Noise 



P5,5 
P9,9 
P5,9 
P9,5 



P6,6 = P7,7 = P8,8 = 8 
PlO.lO = Pll,ll = Pl2,12 
P6,10 = P7,ll = P8,12 = 
Pl0,6 = Pll,7 = Pl2,8 = 



/ 



M34, M4 



P5,5 — P6,6 — Jq^PTJ — f- 
3 

P9,9 = PlO.lO = 16 '^11.11 

3 

P5,9 = P6,10 = 16' ^7,11 = 

3 

P9,5 = Pl0,6 = ig,Pll,7 = 



^,8 



16 



16 



- Pl2,12 — — v^ 
P8,12 = JQ 
Pl2,8 = JQ 



M3^,Mi 
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Pl2,12 — 16 v^ 



_3_ 
16 
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P5,5 = 2'P6,6 = P7,7 = / 
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P5,9 = 2'P6,10 = P7,ll = 
P9,5 = 2;PlO,6 = Pll,7 = 
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= Pl2,12 = 

P8,12 = 
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/ 



M34, M4 
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0,P6,6 — f )P7,7 
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)Pll,ll — Pl2,12 — 



/ 



0,P6, 

0,pm 



10 



6 



2)P7,11 
2'Pll.7 



P8,12 
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M34, M4 



P5,5 
P9,9 

P5,9 
P9,5 



P6,6 = 0,P7,7 
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2?P8,8 
0?Pll,ll = 55Pl2,12 



P6,10 
Pl0,6 



0,P7,11 
0,Pll,7 



2)P8,12 
2'Pl2,8 



/ / S.AR. -^ M34 

S.R.+S.AR. -> M4 








P5,5 
P9,9 

P5,9 
P9,5 



P6,6 — P7,7 — 0, ps,8 — 2 
PlO,10 = Pll,ll = 0,pi2,12 



/ 



P6,10 
Pl0,6 



P7,ll 
Pll,7 



0,P8, 

0,P 



12 



12,8 



/ N.S^M34,M4 

S.R. -^ M234, M2 



Dependence on state preparation of noise induced effects in multiquhit systems 23 

References 

[1] Nielsen M A and Chuang I L 2000 Quantum Computation and Quantum Information (Cambridge: 

Cambridge University Press) 
[2] Schlosshauer M 2007 Decoherence and the Quantum-to- Classical Transition (Berlin; Springer) 
[3] Lidar D A and Whaley K B 2003 arXiv:quant-ph/0301032 
[4] Lidar D A, Chuang I L and Whaley K B 1998 Phys. Rev. Lett. 81 25942597 
[5] Knill E and Lallamme R 1997 Phys. Rev. A 55 900 

[6] Franco R L, Bellomo B, Maniscalco S and Compagno G 'a rXiv:1205. 64191 
[7] Misra B and Sudarshan EGG 1977 18 756 

[8] Maniscalco S, Francica F, Zaffino R L, Gullo N L and Plastina F 2008 Phys. Rev. Lett. 100 090503 
[9] Francica F, Plastina F and Maniscalco S 2010 Phys. Rev. A. 82 052118 
[10] Ando B and Graziani S 2000 Stochastic Resonance: Theory and Applications (Dordrecht: Kluwer 
Academic Publishers) 

Rivas A, Oxtoby N P and Huelga S F 2009 Eur. Phys. J. B 68 51 

Wellens T, Shatokhin V and Buchleitner A 2004 Rep. Prog. Phys. 67 45 

Ghikas D P K and Tzemos A C 2012 Int. J. Quantum Inform. 10 1250023 

Yu T and J. H. Eberly H 2005 Phys. Rev. B 68 165322 

Yu T and J. H. Eberly H 2004 Phys. Rev. Lett. 93 140404 

Yu T and Eberly H 2006 Opt. Commun. 264 393 

Yu T and Eberly H 2007 arXiv:quant-ph/0503089>^3 



[11 

[12; 
[13; 

[14 

[15; 
[16; 

[17 

[18; 

[19 

[20 
[21 
[22; 
[23; 
[24 
[25; 
[26; 
[27 

[28 
[29; 
[30 
[31 
[32; 
[33 
[34; 
[35 
[36 



Yu T and Eberly H 2007 arXiv:quant-ph/0703083 



Yu T and Eberly H 2009 ■arXiv:09lO. 1396' 
Terra Cunha M. O. 2007 New J. Phys. 9 237 

Wang J, Batelaan H, Podany J and Starace A F 2006 J. Phys. B 39 4343 
Jordan P and Wigner E 1928 Z. Phys. 47 631 61 
Lieb E, Schultz T and Mattis D 1961 Ann. Phys. 16 407 61, 67 

Parkinson J B and Farnell D J J 2010 An Introduction to Quantum Spin Systems (Berlin: Springer) 
Rao K R K and Kumar A 2012 Int. J. Quantum Inform. 10 1250039 
Yeo Y, Liu T, Lu Y-E and Yang Q-Z 2005 J. Phys. A: Math. Gen. 38 3235 
Luczka J. 1991 Czech. J. Phys. 41 289-292 

Horodecki R, Horodecki P, Horodecki M and Horodecki K 2005 Rev. Mod. Phys. 81 865942 
Horodecki M, Horodecki P and Horodecki R 1998 Phys. Rev. Lett. 80 5239 
Mintert F, Carvalho A, Kus M and Buchleitner A 2005 Phys. Rep 415 4 207-259 
Hill S and Wooters W K 1997 Phys. Rev. Lett. 78 5022-5025 
Wooters W K 1998 Phys. Rev. Lett. 80 2245-2248 
Rau A R P 2009 J. Phys. A: Math. Theor 42 412002 
Ah M, Rau A R P and Alber G 2010 Phys. Rev. A 81 042105 
Quesada N, Al-Qasimi A and James D F V 2012: arXiy: 1207.3689l 

Breuer H P and Petruccione F 2006 The Theory of Open Quantum Systems (Oxford: Oxford 
University Press) p 160-166 



